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GAMMA-RAY BURSTS MAY BE BIASED TRACERS OF STAR FORMATION 
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ABSTRACT 

Based on a simulation of galaxy formation in the standard cosmological model, we suggest that a 
consistent picture for Gamma-Ray Bursts and star formation may be found that is in broad agreement 
with observations: GRBs preferentially form in low metallicity environments and in galaxies substantially 
less luminous that L» . We find that the computed formation rate of stars with metallicity less than 0.1 Z Q 
agrees remarkably well with the rate evolution of Gamma-Ray Bursts observed by Swift from z — to 
z = 4, whereas the evolution of total star formation rate is weaker by a factor of about 4. Given this 
finding, we caution that any inference of star formation rate based on observed GRB rate may require 
a more involved exercise than a simple proportionality. 

Subject headings: stars: abundances — supernovae: general — galaxies: formation — cosmology: 
theory 



1. INTRODUCTION 

The intriguing observational linkage between long dura- 
tion Gamma-Ray Bursts (GRBs) and core-collapse super- 
novae (e.g., Stanek et al. 2003; Hjorth et al. 2003) suggests 
that the progenitors of GRBs may be very massive stars. 
This possible connection was predated by a proposed uni- 
fied picture (Cen 1998). As such, it has been hoped that 
GRBs may be a good tracer of cosmic star formation (Wi- 
jers et al. 1998; Totani 1999; Lamb & Reichart 2000; Blain 
& Natarajan 2000; Porciani & Madau 2001; Daigne et al. 
2006; Coward 2006; Le & Dermer 2007; Li 2007). How- 
ever, recent observations indicate that typical GRBs may 
prefer relatively low metallicity environments (~ 0.1 Zq) 
(Fynbo et al. 2003; Le Floch et al. 2003; Christensen et 
al. 2004; Fruchter et al. 2006; Stanek et al. 2006) and 
host galaxies significantly less luminous than L. t (Fruchter 
et al. 1999,2006), although there is evidence that the ac- 
tual spread in metallicity may be wide (Berger et al. 2006; 
Prochaska 2006; Wolf & Podsiadlowski 2007). 

The aim of this Letter is to first address the issue of 
consistency of GRB environment with respect to metal- 
licity and galaxy luminosity, i.e., the galaxy luminosity- 
metallicity relation, in the context of detailed simulation 
of galaxy formation in the standard cosmological model. 
Then, we make predictions on the evolution of GRB rate 
with redshift and highlight a possible dramatic difference 
between overall star formation history and GRB rate his- 
tory, if GRBs are not an unbiased tracer of star forma- 
tion. In particular, if GRBs are predominantly produced 
by stars with metallicity < 0.1 Z , the GRB rate is ex- 
pected in our model to rise obstinately from z = to 
z ~ 5 by a factor of ~ 100, when it flattens out towards 
higher redshift, whereas the overall star formation rate 
rises rapidly only from z — to z ~ 3 and is roughly flat 
from z > 3 until z ~ 7. The evolution of GRB rate with 
redshift is thus expected to be stronger than that of star 
formation. 

2. EVOLUTION OF GRB AND STAR FORMATION RATES 



Observations seem to indicate that GRB galaxy hosts 
are preferentially dwarf galaxies of about 0.1L* at low red- 
shift (Fruchter et al. 1999,2006). This would be surprising, 
if GRB rate is directly proportional to total star formation 
rate, because the latter is known to peak in significantly 
larger galaxies (see, for example, Figure 1 below). The im- 
plication is that GRB rate is not exactly proportional to 
the overall star formation rate. On the other hand, anal- 
ysis of the metallicity of GRB progenitors suggests that 
GRB progenitors tend to have relatively lower metallicity 
of ~ 0.1 Z than that of typical forming stars at low red- 
shift (Fynbo et al. 2003; Le Floch et al. 2003; Christensen 
et al. 2004; Fruchter et al. 2006; Stanek et al. 2006). We 
would like to ask the following question: are these two ob- 
servational facts consistent in the galaxy formation model 
in the standard cosmological model? 

Figure 1 shows the distribution of SDSS U-band light as 
a function of the stellar mass of galaxies, with the galaxies 
being divided into two subgroups according to the mean 
metallicity of galaxies. The results are based on a cold 
dark matter cosmological simulation with galaxy forma- 
tion in a cold dark matter universe (Nagamine et al. 2006; 
Cen & Ostriker 2006; Cen & Fang 2006) with the following 
essential parameters: SIm = 0.31, = 0.048, ft a = 0.69, 
cr 8 = 0.89, H = 100/ikm s" 1 Mpc" 1 = 69 km s" 1 Mpc" 1 
and n s = 0.97. The simulation box size is 85ft. _1 Mpc co- 
moving with a number of cells of 1024 3 , giving a cell size of 
83 h _1 kpc comoving and dark matter particle mass equal 
to 3.9 x 10 8 /i _1 Mq. Given a lower bound of the tempera- 
ture for almost all the gas in the simulation of T ~ 10 4 K, 
the Jeans mass is ~ 10 10 M for mean-density gas, which 
is comfortably larger than our mass resolution. Galaxies 
are produced using a grouping scheme HOP (Eisenstein 
& Hutl998) (see Nagamine et al. 2001 for details), which 
gives a catalog of galaxies each with stellar mass, mean 
stellar metallicity, luminosities in all SDSS bands, etc. 

We see a general trend in Figure 1 that larger galaxies 
tend to have higher metallicity, albeit a significant disper- 
sion exists (not shown here). Specifically, it is seen that, 
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if the stellar metallicity is, say, typically lower than 0.1 Zq 
for GRBs progenitors, galaxies of stellar mass 10 10 M is 
expected to make the dominant contribution, whereas the 
overall star formation rate is peaked in galaxies of stel- 
lar mass 10 11 — 1O 12 M0. Apparently, this is in broad 
agreement with observations of GRB being concentrated 
in dwarf galaxies of metallicity < 0.1 Zq. This is reassur- 
ing in that our simulation that has been shown to produce 
consistent results compared to the real universe in many 
other aspects appears to be in broad agreement with obser- 
vations with regard to the general luminosity-metallicity 
trend (e.g., Kobulnicky & Kewley 2004) and with the 
observed preference of GRB hosts of typical metallicity 
~ 0.1 Zq and typical luminosity < 0.1L*. 
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Fig. 2. — shows histories of three rates: total star formation (dot- 
ted curve), star formation with metallicity less than 0.3 Zq (dashed 
curve) and less than 0.1 Zq (solid curve), respectively. 

We note that, while the overall star formation seen in 
Figure 2 is higher by a factor of ~ 30 at z — 4 than at 
2 = 0, the ratio of GRB rate at z = 4 to that at z = is 
~ 120, if GRB progenitors predominantly have metallicity 
less than 0.1 Z , roughly a factor of four larger. This result 
is, curiously, in remarkable agreement with recent obser- 
vations of Kistler et al. (2007) where they suggest that the 
GRBs observed by Swift has an enhanced evolution by a 
factor of ~ 4 from z — to z = 4 compared to the overall 
star formation rate. 



Fig. 1. — shows the distribution of U band light as a function 
of the stellar mass of galaxies, for two groups of galaxies, one with 
metallicity greater than 10% of solar metallicity (dashed histogram) 
and the other with metallicity lower than 10% of solar metallicity 
(solid histogram). 

Having verified that our simulation is able to reproduced 
the general luminosity-metallicity relation for galaxies, we 
will take a step further to infer the rate evolution for 
GRBs, given the computed star formation rate and stellar 
metallicity as a function of redshift. Since our simulated 
galaxies are composed of thousands to millions of "stel- 
lar" particles, which typically resemble globular clusters, 
we will use the "resolved" metallicity of individual "stel- 
lar" particles. Usually, there is a wide dispersion in stellar 
metallicity among the "stellar" particles within an individ- 
ual galaxy, reflecting complicated star formation history of 
each galaxy. 

Figure 2 shows histories of three rates: total star forma- 
tion (dotted curve), star formation with metallicity less 
than 0.3 Z Q (dashed curve) and less than 0.1 Z Q (solid 
curve), respectively. We see that, if GRBs are preferen- 
tially hosted by stellar progenitors with metallicity lower 
than 0.1 Zq, one should expect to see their rate to rise 
approximately exponentially as a function of redshift from 
z — until z ~ 5. This is contrasted with a flattening 
of the overall star formation rate at a much lower redshift 
z ~ 2 accompanied by only a modest rise (a factor of 2) 
from z ~ 2 to z ~ 5. 
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Fig. 3. — shows the cumulative distribution of the observed (steps) 
and predicted (lines) GRB rates, normalized to the total number of 
bursts between z=0 and z=4. See Fig. 2 for line styles. Data 
shows the distribution of the 36 Swift GRBs taken from Butler et 
al. (2007). The star formation rate of all-galaxy sample was nor- 
malized to the analytic values of Hopkins & Beacom (2006), and 
the star formation rates for the other two samples was normalized 
by the same factor. 

In the spirit of making as a direct comparison as possi- 
ble with observations, we shows the cumulative distribu- 
tion of the observed (steps) and predicted (curves) GRB 
rates, normalized to the total number of bursts between 
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z = and z = 4 in Figure 3, following Kistler et al. (2007). 
Here, in order to make a direct comparison with Kistler et 
al. (2007) and to remove uncertainties in the overall star 
formation rate in our simulation, we normalize the star 
formation rate of the all-galaxy sample in our simulation 
to the analytic values of Hopkins & Beacom (2006), as 
did Kistler et al. (2007). Then, the star formation rates 
for the other two (lower metallicity) samples are adjusted 
multiplicatively by the same factor. It is seen that, if one 
simply assumes that the total GRB rate is proportional 
to the star formation rate of metallicity less than 0.1 Z , 
the observed evolution of GRB rate from z = to z = 4 
is reproduced to a high degree, clearly visible in the good 
agreement between the solid curve and the step curve in 
Figure 3. 

3. DISCUSSION 

Taken all the observational facts together along with our 
theoretical results, a broadly consistent picture appears to 
emerge: GRBs preferentially form in low metallicity envi- 
ronments and (because of the luminosity-metallicity rela- 
tion) in galaxies substantially less luminous that L*, and 
GRB rate evolves more strongly than the overall star for- 
mation rate. Evidently, the observed evolution of GRB 
rate from z = to z = 4 can be explained, if GRBs are 
primarily produced by massive stars with metallicity less 
than 0.1 Zq. Including higher metallicity stars would pro- 
duce GRB rate evolution from z = to z — 4 that is 
inconsistent with observations. Nevertheless, this overall 
picture would also be consistent with a theoretical prefer- 
ence or possibly requirement of low metallicity for GRB 
progenitors in the context of "collapsar" models (MacFay- 
den & Woosley 1999; Woosley & Heger 2006). 

What is implicitly assumed is that the stellar initial 
mass function (IMF) has remained the same over the red- 
shift range considered. In other words, whatever metallic- 
ity dependence of GRB rate may have, this dependence is 
assumed not to evolve with redshift. One should note that 
it is not fully known observationally or understood theo- 
retically how the IMF evolves with time. Therefore, addi- 
tional possible evolutionary effect of IMF would add an- 
other layer of complexity to this issue. It is often thought 
that lower metallicity environment might favor formation 
of more massive stars. If GRB progenitors are massive 
stars, this would then translate to the expectation that ad- 
ditional effect due to an evolving IMF may further steepen 



the evolution of the GRB rate with redshift. This, how- 
ever, is not required or borne out in our analysis. Our 
results thus imply that the evolution of IMF from z = 4 
to z = 0, if any, appears to be weak. 

If we place this result in a larger context of star forma- 
tion over the entire cosmic history, one might come to the 
conclusion that, while there may be a dramatic transition 
of IMF from Population III metal-free stars (Nakamura& 
Umemura, M. 2002; Abel, Bryan & Norman 2002; Bromm, 
Coppi, & Larson 2002) to Population II stars at some high 
redshift (e.g., Cen 2003; Trac & Cen 2007), further evolu- 
tion of IMF at lower redshift may be modest, in the sense 
that the mass fraction of high mass stars that are presum- 
ably GRB progenitors of the total stellar mass remains 
relatively constant. 

4. CONCLUSIONS 

We utilize a simulation of galaxy formation in the stan- 
dard cosmological model that has been shown to produce 
results consistent with extant observations of galaxy for- 
mation (e.g., Nagamine et al. 2006) to shed light on the 
relation between GRB rate and star formation rate. We 
find that a consistent picture for Gamma-Ray Bursts and 
star formation that is in broad agreement with observa- 
tions would emerge, if GRBs preferentially form in low 
metallicity environments and in galaxies substantially less 
luminous that L*. Because of the increase of metallicity 
with cosmic time, GRB rate consequently evolves more 
strongly with redshift than the overall star formation rate. 
We find that the observed evolution of GRB rate from 
z = to z = 4 can be explained, if GRBs are primarily 
produced by massive stars with metallicity less than 0.1 Z Q , 
whereas an inclusion of stars with metallicity as high as 
0.3 Z© yields GRB rate evolution from z = to z = 4 
inconsistent with observations. 

Therefore, we reach the conclusion that GRBs may not 
be a good tracer of cosmic star formation, especially over 
a long timeline. As a result, a simple inference of star for- 
mation rate or its derived quantities such as the ionizing 
photon production rate at high redshifts, based on the ob- 
served GRB rate, should be done with caution and may 
require careful calibrations. 

We thank Ken Nagamine for providing simulated galaxy 
catalogs. We gratefully acknowledge financial support by 
grants AST-0407176 and NNG06GI09G. 
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